%1*2^3+1*2^2+1*2^1+1*2^0+1*2^-1+1*2^-2+1*2^-3+1*2^-4
clc
clear

r = 2;
L = 4;
M = 4;
a = -1;
b = 1;
system_no = 1;
system1(system_no,r,L,M,a,b);

r = 4;
L = 4;
M = 4;
a = -2;
b = 2;
system_no = 2;
system1(system_no,r,L,M,a,b);

r = 4;
L = 4;
M = 4;
a = -3;
b = 3;
system_no = 3;
system1(system_no,r,L,M,a,b);

r = 7;
L = 3;
M = 3;
a = -4;
b = 6;
system_no = 4;
system1(system_no,r,L,M,a,b);